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ABSTRACT 

Using data from the Fermi Large Area Telescope (LAT), we report the first clear 7 -ray measurement 
of a delay between flares from the gravitationally lensed images of a blazar. The delay was detected 
in B0218+357, a known double-image lensed system, during a period of enhanced 7 -ray activity with 
peak fluxes consistently observed to reach > 20 — 50 x its previous average flux. An auto-correlation 
function analysis identified a delay in the 7 -ray data of 11.46 ±0.16 days (lcr) that is ~ 1 day greater 
than previous radio measurements. Considering that it is beyond the capabilities of the LAT to 
spatially resolve the two images, we nevertheless decomposed individual sequences of superposing 
7 -ray flares/delayed emissions. In three such ~ 8 — 10 day-long sequences within a ~ 4-month span, 
considering confusion due to overlapping flaring emission and flux measurement uncertainties, we 
found flux ratios consistent with ~ 1 , thus systematically smaller than those from radio observations. 

During the first, best-defined flare, the delayed emission was detailed with a Fermi pointing, and we 
observed flux doubling timescales of ~ 3 — 6 hrs implying as well extremely compact 7 -ray emitting 
regions. 

Subject headings: Galaxies: active — gravitational lensing: strong — gamma rays: galaxies — quasars: 
individual (B0218+357) 
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1. INTRODUCTION 

B0218+357 was discovered with the NRAO 140- 
ft telescope in its strong source survey (S3 0218+35; 
Pauliny-Toth & Kellermann 1972). Later radio imag- 
ing revealed it to be a gravitationally lensed blazar with 
the smallest separation double-image known (335 milli- 
arcseconds), and an Einstein ring with a similar an- 
gular diameter (O’Dea et al. 1992; Patnaik et al. 1993). 
The lens galaxy is at redshift z = 0.6847 (Browne et al. 
1993), and the blazar was later securely measured at 
2 = 0.944 ± 0.002 (Cohen et al. 2003). 

Shortly after the lens discovery, Corbett et al. (1996) 
measured a time delay (Refsdal 1964) A t r = 12 + 3 days 
(lcr quoted throughout unless otherwise specified) at ra- 
dio wavelengths, using the VLA to spatially separate and 
monitor the polarization variability in its leading brighter 
A (western) and fainter B (eastern) images. Later in- 
dependent (but contemporaneous) dual-frequency VLA 
observations further refined the delay, Af r = 10.5 ± 0.2 
(Biggs et al. 1999) and 10.1 ± 0.8 days (Cohen et al. 
2000). Interestingly, Eulaers & Magain (2011) analyzed 
the latter’s measurements and found two possible delays, 
Af r = 9.9to:g or 11-8 ± 2.3 days. Although these delays 
span a narrow range, Af r ~ 10 — 12 days, because of the 
differing assumptions and analysis techniques employed 
in these works, there remains some debate on how to best 
derive their uncertainties. 

B0218+357 is also a y-ray source detected by the 
Fermi Large Area Telescope (LAT; Atwood et al. 2009) 
with an average flux 30 + 7 = (1.00 ± 0.07) x 10 -7 pho- 
tons cm -2 s _1 over its first 2-years of observations 
(2FGL J0221. 0+3555; Nolan et al. 2012). Its steep spec- 
trum at > 100 MeV energies (photon index, T = 2.28 ± 
0.04) and overall spectral energy distribution are typical 
of an otherwise normal y-ray emitting flat-spectrum ra- 
dio quasar (e.g., Ghisellini et al. 2010). While y-ray data 
lack the necessary spatial resolution to separate lensed 
images, such blazars display their most dramatic vari- 
ability in y-rays, and the LAT’s all-sky monitoring could 
give it a distinct advantage over lower-frequency imaging 
observations in parameterizing lensed systems. Indeed, 
Atwood (2007) proposed prior to Fermi’s launch that 
the LAT could detect delayed emission from such gravi- 
tationally lensed blazars using integrated lightcurves for 
sufficiently bright y-ray flares. B0218+357 was found to 
be variable in the early LAT observations, though only 
modestly so (Abdo et al. 2010, see Figure 1). 

Bright y-ray flaring from B0218+357 was detected 
with the LAT beginning late 2012 August (Ciprini 2012), 
and a delayed flare was tentatively identified ~ 10 
days later (Giroletti et al. 2012), consistent with the 
radio delay measurements. The blazar then displayed 
even brighter, more sustained flaring activity begin- 
ning September 14, thus prompting a Fermi Target-of- 
Opportunity (ToO) pointed observation (Cheung et al. 
2012) that traced the anticipated delayed emission in 
detail. Two main additional flaring events were subse- 
quently observed in as many months (see Figure 1 for an 
overview). We discuss the y-ray temporal and spectral 
properties of B0218+357 together with the derived time 
lag, flare timescales, and observed flux ratios of the A/B 

30 LAT 7 -ray fluxes are reported at E > 100 MeV throughout. 


images. 

2. LAT OBSERVATIONS AND ANALYSIS 

The Fermi- LAT operates in a default sky-survey mode, 
and over every two ~1.6-hr spacecraft orbits, provides 
observations covering the entire sky. We used LAT ob- 
servations with the P7S0URCE_V6 instrument response 
functions, selecting 100 MeV - 100 GeV events with 
a region of interest (ROI) of radius = 15° centered 
at the B0218+357 radio position, R.A. = 35°. 27279, 
Deck = 35°. 93715 (J2000; Patnaik et al. 1992). The 
maximum zenith angle of 100° was set to minimize the 
contamination from Earth limb photons as well as the 
appropriate gtmktime filter (#3) following the FSSC 
recommendations 31 for the combination of sky-survey 
and pointed observations. The gtlike likelihood in 
the Fermi Science tools (version v9r27pl) was used 
for the spectral analysis, assuming throughout a single 
power-law model for B0218+357 over the selected en- 
ergy range (as in the 2FGL catalog). The background 
model included all 2FGL sources within the ROI as well 
as the Galactic (gal_2yearp7v6_v0.fits) and isotropic 
(iso_p7v6source.txt) diffuse components. 

In generating each lightcurve, the isotropic normaliza- 
tion was left free to vary in each time-bin while the two 
known variable 2FGL sources within a 5° ROI and the 
Galactic normalization were initially fitted over each full 
interval, then fixed at the average fitted values in the 
shorter time-bins. As a convenient reference point, we 
define T = MJD — 56100 days (i.e., T = 0 was 2012 June 
22), the time when y-ray flaring became obvious. Inte- 
grating 1417 days (~ 3.9 years) of LAT observations prior 
to this date gave an average F 7 = (0.83 ± 0.05) x 10~ 7 
photons cm' 2 s _1 , with T = 2.30 + 0.03, consistent with 
the 2FGL value. For context, we generated a 1-week 
binned lightcurve for 5-years of data (2008 August 5 - 
2013 August 6; Figure 1, top) assuming a fixed T = 2.3. 
Besides the modest source activity in early 2009 and 
2010, the pronounced flaring beginning in mid-2012 last- 
ing for ~ 200 days is apparent; thereafter, the source 
quieted again to earlier levels. 

In order to study the flaring activity in detail, we de- 
fined a 265-day interval starting at T = 0 days and gener- 
ated 1-day and 6-hr binned lightcurves. The Fermi ToO 
observations also allowed us to produce a ~ 1.6- hr orbit- 
by-orbit binned lightcurve for the sub-interval covering 
the first delayed flare from 2012 September 24 - Octo- 
ber 1 (T = 94 — 101 days). To search for any possible 
spectral changes, we initially computed the 1-day binned 
lightcurve with the photon index free in the fit. For the 
108 points with the greatest significance (test statistic 32 , 
TS > 25), we found all but four points within 2cr of the 
weighted average value of 2.31 ± 0.02, which in turn is 
consistent with the 3.9-yr average. We thus regenerated 
the 1-day (Figure 1, middle), the 1.6-hr orbit (Figure 1, 
bottom), and the 6-hr binned lightcurves (Figure 2) with 
T = 2.3 fixed. 

3. RESULTS 

31 http: / /f ermi . gsfc.nasa.gov/ ssc/data/analysis/ 
documentation/ Cicerone/Cicerone_Likelihood/ Exposure . 

html 

32 The source significance is equivalent to ~ y/TS , assuming one 
degree of freedom (Mattox et al. 1996). 
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Fig. 1. — LAT lightcurves over wide dynamic range: 1-week bins over the first 5-years of the Fermi mission (top), 1-day bins for a 265-day 
flaring interval (middle), and 1.6-hr orbit bins during the 7-day Fermi ToO (bottom). The pre, post, and three active episodes outlined 
in the middle panel are further detailed in Figure 2. Throughout, flux points (plotted with 1 <t errors when TS > 4 in the bin) and arrows 
indicating 2a upper limits (when TS < 4) are connected by dotted lines. Horizontal dashed lines indicate the 3.9-year average flux prior 
to the flaring interval (top) and the baseline flux during the flaring interval, F 7 = 0.3 X 1 0 fl photons cm -2 s — 1 (middle, bottom). 


3.1. Time Lag 

The B0218+357 7-ray lightcurve appears quite com- 
plex with many peaks and valleys over the ~ 4 months 
from T ~ 60 — 180 days (Figure 2) when the source 
was most active. To search for a time lag, we com- 
puted the auto-correlation function (ACF) for the 6- hr 
binned lightcurve up to lag values of half of the total 
defined 265-day flaring interval. This evenly sampled 
lightcurve consisted of 1057 measurements with three 
missing data points due to exposure gaps. The ACF 
was therefore computed both by a standard (after inter- 
polating the three missing points) and a discrete routine 
(Edelson & Krolik 1988). The two procedures gave al- 


most identical results and the ACF is shown in Figure 3. 
A single prominent correlation peak is apparent between 
the time lag range of 11 — 12 days. The peak’s signifi- 
cance is 9cr with respect to the measurement noise and 
comparing it to the height above the ACF “background.” 
Fitting a Gaussian function to this peak, we estimated a 
best-fit value, A t 7 = 11.46 ± 0.16 days (lcr). Uncertain- 
ties were estimated by a model independent Monte Carlo 
method (Peterson et al. 1998) accounting for the effects 
of measurement noise and data sampling. The time lag 
does not match any known period observed with the 
LAT (Ackermann et al. 2012; Corbet et al. 2012). Be- 
cause the 7-ray flaring was so pronounced especially from 
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Fig. 2. — LAT lightcurve in 6-hr bins from 2012 July 22 - December 24 detailing the pre, post, and three main episodes 
(cf., Figure 1, middle), sub-divided into ~ 8 — 10 day-long flares and corresponding delays (asterisks mark outlying sharp 
features; see text). Each panel spans 55-days, with adjacent panels overlapping by 5-days on each side. Horizontal dashed lines 
indicate the baseline flux during the flaring interval (Figure 1). 


T ~ 84 — 155 days, and appears to be broadly divided 
into three ~ 8 — 10 day long flare/delay sequences (Sec- 
tion 3.2), this could induce other smaller enhancements 
in the ACF over the studied interval. 

As a cross-check of the lag derived from the full flar- 
ing interval 7-ray data, discrete ACFs were computed 
for two segments from T = 0 — 110 and T = 110 — 265 
days. The lags obtained from Gaussian fits to the peaks 
were, A t 7 = 11.52 ± 0.31 and Af 7 = 11.38 ± 0.28 days, 
respectively, confirming the delay value and small uncer- 
tainty for the full interval, thus indicating that we ob- 
tained a robust measurement with the LAT. The small 
uncertainty in Af 7 is competitive with the best deter- 
mined radio measurements for B0218+357 although the 
former is marginally larger by At 7 — Af r = 1.0 ± 0.3 
and 1.4 ± 0.8 days (lu) than the Biggs et al. (1999) and 


Cohen et al. (2000) values, respectively, but consistent 
with the Eulaers & Magain (2011) values. If the radio/7- 
ray delays are intrinsically different due to an offset be- 
tween the respective emitting regions, the implied offset 
in a singular isothermal sphere lens model is ~ 70 pc 
(projected) for a ~ 10% difference in the time delay. 
This seems extreme considering such offsets are on av- 
erage ~ 7 pc in other blazar jets (e.g., Pushkarev et al. 
2010), and may rather suggest the uncertainty in the ra- 
dio delay was underestimated (Section 1). 

3.2. Flare Timescales 

Utilizing the 7-ray delay measurement, we can broadly 
identify three sets of flare/delay episodes in the 6-hr 
binned LAT lightcurve of B0218+357 (Figure 2). The 
pre-flare times were what triggered the initial excitement 
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Fig. 3. — Auto-correlation function computed for the 6-hr binned 
LAT lightcurve of the 265-day flaring interval. The inset zooms in 
around the best-fit indicated lag peak. 

in late 2012 August and are now detailed as a 6-hr flare 
at T ~ 50 days (with a corresponding delayed signal 
11.5 days later) and a doublet of 6 — 12-hr flares 1 day 
apart beginning at T ~ 65 day. In the doublet, only 
the first flare showed a clear delayed flare 11.5 days later 
while the second shows no similar corresponding delayed 
(or 11.5-day prior) feature; microlensing (see below) or 
a relatively large variation in the magnification ratio are 
possible explanations. 

The first bright 7-ray sequence began at T = 84 day 
with the best-defined flaring structure with observed 
fluxes, ~ (2 — 5) x 10 -6 photons cm -2 s _1 over eight 
consecutive 6-hr bins, followed by a sharp drop and sub- 
sequent rise in one day. The Fermi ToO observation be- 
gan 10 days later and the anticipated delayed emission 
mirrored the initial flare with the rise and peak sepa- 
rated by 1 day and all features well-matched 11.5 days 
later. We broadly identified two subsequent (2nd and 
3rd) ~ 8 — 10 day duration 7-ray flaring sequences, but 
these were more difficult to disentangle because of super- 
posing flares in the integrated lightcurves. The post-flare 
intervals showed lower fluxes, comparable to the pre-flare 
emission states. 

In Figure 2, the observed variability timescales (dou- 
bling and halving), t var , during the first and subsequent 
two flaring episodes are securely less than the 6-hr bin- 
ning. Doubling timescales as short as 2-orbits (~ 3 
hrs) are further suggested in the orbit-by-orbit binned 
lightcurve from the Fermi ToO pointing of the first 
delayed flare (Figure 1, bottom). Such timescales are 
amongst the fastest well-constrained 7-ray variability in 
a blazar observed with the LAT (Tavecchio et al. 2010; 
Abdo et al. 2011) and constrain the 7-ray emission re- 
gion diameter, d < 2ct vav / (l + z)<6x 10 14 cm, mod- 
ulo the unknown Doppler beaming factor. Assuming an 
h = H 0 / {100 kms" 1 Mpc" 1 ) = 0.71 (fi M = 0.27, fl A = 
0.73) cosmology, this translates to an angular diameter 
« 30 nano-arcseconds, ~ 10 4 x smaller than the best ra- 
dio size constraint (Mittal et al. 2007). Microlensing is 
thus an important factor in interpreting our 7-ray re- 
sults because the smaller the structures, the larger the 
expected variability of magnification. 

3.3. Flux and Magnification Ratios 


Adopting the 7-ray delay, we compared the 6-hr binned 
lightcurves for the three main flaring episodes with the 
observations shifted by —11.46 days and computed the 
observed ratios between corresponding flux pairs, retain- 
ing only ratio values > 2x their uncertainties (Figure 4). 
The first sequence appears to show the clearest corre- 
spondence between features in the two lightcurves, with 
only minor deviations about the weighted average flux ra- 
tio 1.3 ± 0.1. By subtracting a baseline, F 1 = 0.3 x 10~ 6 
photons cm -2 s _1 (the minimum observed flux during 
the overall flaring interval), we can further estimate a 
corresponding magnification ratio in 7-rays s=s 1.3, con- 
sistent with the flux ratio. The average ratios for this 
first sequence seem to imply the brighter A image led 
the B image in 7-rays, as observed in the radio. More 
conservatively however, given the large uncertainties in 
the individual measurements, the flux ratios appear con- 
sistent with unity. Moreover, for the subsequent 2nd and 
3rd sequences, the correspondences between the flare and 
delayed emissions were less clear. Sharp and more scat- 
tered changes in the paired flux ratios were apparent, 
including values < 1 (which would imply a fainter lead- 
ing A image). We interpret this as an artifact due to 
contamination from superposing flares after the source 
has already entered a very active phase. This confusion 
in the integrated lightcurves prevents us from reliably 
determining magnification ratios, and how variable this 
quantity may have been. 

The flux ratio measured in 7-rays is smaller than in 
the radio. Biggs et al. (1999) found a small, but statisti- 
cally significant frequency dependence in the flux ratios, 
3.57 ± 0.01 (8 GHz) and 3.73 ± 0.01 (15 GHz), while 
Cohen et al. (2000) found similar values but with larger 
uncertainty, 3. 2+^4 (8 GHz) and 4.3to's (15 GHz). Fre- 
quency dependence in the flux ratios of the two radio im- 
ages and their observed substructures could be possibly 
due to free-free absorption and scattering from a molecu- 
lar cloud in the lens galaxy (Mittal et al. 2007). We note 
that the radio and 7-ray observations are not simulta- 
neous and magnification ratios could be variable with 
time. Further complicating such comparisons are open 
questions in blazar jet studies, i.e., the radio and 7-ray 
emitting regions need not coincide, with the latter likely 
more compact (Section 3.2), and whether successive 7- 
ray flares originated in a single emission zone or from 
separate relativistically moving dissipation regions along 
the jet. Excursions could also be due to intrinsic changes 
in the magnification ratios or microlensing from the rela- 
tive motion of the source seen through a clumpy lensing 
galaxy. Indeed, microlensing in the context of extremely 
compact 7-ray emission zones (Torres et al. 2003) could 
explain the single 6-hr flare points that do not have cor- 
responding lags (marked with asterisks in Figure 2), al- 
though fast superposed flares are also a possibility. Note 
that in optical and infrared observations, the B image ap- 
pears brighter than the A image, i.e., reversed from the 
radio situation, and this is likely due to a combination of 
extinction of the A image and microlensing (Falco et al. 
1999; Jackson et al. 2000). 

4. DISCUSSION AND CONCLUSIONS 

Our detection of a gravitational lens time delay, Af 7 = 
11.46 ± 0.16 days, in the LAT observations of blazar 
B0218+357 has some interesting potential implications 


6 


Cheung et al. 



Fig. 4. — Top panels show the 6-hr binned lightcurves (Figure 2) for the three flares (filled blue) and delayed emission shifted by —11.46 
days (open red). Bottom panels show individual observed flux ratios (dashed line drawn at ratio = 1 for reference) in the corresponding 
upper panels; error bars are symmetric and the third panel was cropped in order to display a common range. 


for future 7-ray studies. Foremost, the LAT detection 
of a 7-ray gravitational lens flaring event in B0218+357 
suggests that such a measurement is possible in other 
blazars. In particular, gravitational lenses found in sur- 
veys of flat-spectrum radio sources (Browne et al. 2003; 
Winn et al. 2000) comprise a relevant sample as these 
form the basis of candidate 7-ray blazar catalogs (e.g., 
Healey et al. 2007). There are ~ 20 gravitational lenses 
from these surveys out of > 10 4 radio sources stud- 
ied with £ 30 mJy at 8 GHz and so far, the two ra- 
dio brightest are detected 7-ray sources, PKS1830-211 
(below) and B0218+357 (out of ~ 10 3 known 7-ray 
blazars; Nolan et al. 2012). The other fainter lensed 
systems are typically less variable at radio frequencies, 
making delay measurements difficult (e.g., Jackson 2007; 
Eulaers & Magain 2011) and while they are not yet re- 
ported 7-ray sources, the all-sky monitoring of Fermi- 
LAT will allow the detection of short-timescale flaring 
7-ray activity in which to attempt delay measurements. 
Importantly, 7-ray measurements constrain lens parame- 
ters free of propagation effects like scintillation (Heeschen 
1984; Lovell et al. 2008) that can hamper radio delay at- 
tempts (Winn et al. 2004), although microlensing may 
be an important limiting factor because 7-ray emitting 
regions are expected to be more compact than in the 
radio. 

The case of B0218+357 appears to be the first clear 
case of a 7-ray detected gravitational lens time delay for 
any astrophysical system. Previously, 7-ray flaring from 
the gravitationally lensed z = 2.507 blazar PKS1830-211 
was detected with the Fermi - LAT (Ciprini 2010) with a 
claimed delay, A t 1 = 27.1 ± 0.6 days (Barnacka et al. 

2011) consistent with the radio measurement, At r = 
26^5 days (Lovell et al. 1998). Subsequent analysis of 
more LAT data, including several prominent flares, did 
not confirm the 7-ray delay (Abdo et al. 2013). If the 7- 
ray delay in PKS1830-211 is assumed to be the same as 


the radio-measured delay, the non-detection of delayed 
7-ray flares implies a magnification ratio in 7-rays much 
larger ( ^ 6) than that observed in the radio (1.52 ±0.05; 
Lovell et al. 1998), thus opposite of what we observed in 
B0218±357. With only two examples studied, no trend 
is clear. However, if microlensing effects can be disentan- 
gled (and in fact utilized as additional constraints on the 
emitting region source size) , magnification ratios in radio 
and 7-ray arising from spatially distinct emission regions 
may be utilized as a probe of differing multi-frequency 
jet structures (see Martf-Vidal et al. 2013). 

A time delay due to gravitational lensing of a back- 
ground source by a foreground object can constrain 
Hubble’s parameter (Refsdal 1964). The original lens 
model for B0218±357 (Biggs et al. 1999) predicted a de- 
lay, At = 7.2+2S)h~ l days (95% confidence). Utiliz- 
ing improved localization of the lensing galaxy, the de- 
lay model uncertainty was reduced to 6.0% (York et al. 
2005), thus deriving h = 0.70 ± 0.05, assuming the often 
quoted Biggs et al. (1999) measured radio delay (cf., Sec- 
tion 3.1). Adopting the York model for our independent 
7-ray measured delay results in h = 0.64 ± 0.04, where 
this quoted uncertainty is due only to the time delay es- 
timate and the statistical uncertainty in the mass model. 
Systematic errors in the modeling, and additional un- 
certainty due to line-of-sight structures (e.g., Suyu et al. 

2012) will likely significantly increase this. Nevertheless, 
it is interesting that the LAT time delay brings the esti- 
mated value of Hubble’s constant down, towards the low 
end of modern measurements (e.g., Planck Collaboration 

2013) . An underdense environment would require this 
inferred h value to increase; including external lensing 
effects in future cosmographic analyses might be impor- 
tant in this system. Moreover, since the radio and 7-ray 
emission regions are likely not co-spatial, the assumed 
radio-derived time-delay function values may be inaccu- 
rate. A fully self-consistent joint modeling of the radio 
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and 7-ray source is needed to resolve this. If the LAT 
can measure a lag in the 7-ray lightcurve of one of the 
previously known systems with wider separation or in a 
new example (below), this can give independent 7-ray 
based constraints on Hubble’s constant. 

One exciting result would be the detection of a lens 
delay in a flaring 7-ray source that is not yet identified 
as a gravitationally lensed system at radio wavelengths 
or otherwise. These could possibly be lensed image pairs 
with flat-spectrum radio sources at smaller separations 
than in the 0.2" resolution of VLA surveys (references 
above). Similar radio lens surveys in the southern hemi- 
sphere are not yet as complete (e.g., Prouton et al. 1999), 
so a 7-ray delay signature in their LAT lightcurves could 
betray the presence of a previously unknown lens sys- 
tem. Such a strategy has been proposed for future wide- 
held optical surveys (Pindor 2005), and the discovery 
potential of the LAT in 7-rays should now be recog- 
nized. Furthermore, with the different flux ratios at ra- 
dio and 7-ray wavelengths, and possible variability of 
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